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Abstract 

The transverse momentum distribution of low-mass Drell-Yan pairs is cal- 
, culated in QCD perturbation theory with all-order resummation. We argue 

that at LHC energies the results should be reliable for the entire transverse 
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momentum range. We demonstrate that the transverse momentum distribu- 
tion of low-mass Drell-Yan pairs is an advantageous source of constraints on 



Qh! the gluon distribution and its nuclear dependence. 



I. INTRODUCTION 

Dilepton production in hadronic collisions is an excellent laboratory for the investiga- 
tions of strong interaction dynamics. This channel provides an opportunity for discovery 
of quarkonium states and a clean process for the study of parton distribution functions 
(PDF). In the Drell-Yan process, the massive lepton-pair is produced via the decay of an 
intermediate Z° boson or a virtual photon 7* with mass M. When M ~ Mz, the high mass 
dilepton production in heavy ion collisions at LHC energies is dominated by the Z° channel 
and is an excellent hard probe of QCD dynamics [1]. In this letter, we demonstrate that the 
transverse momentum distribution of low-mass (Aq CD <C M <C M z ) dilepton production at 
LHC energies is a reliable probe of both hard and semihard physics at LHC energies and 
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is an advantageous source of constraints on gluon distribution in the proton and in nuclei. 
In addition, it provides an important contribution to dilepton spectra at the LHC, which is 
the appropriate channel to study J/ip, heavy quarks etc. 

II. FULL P T SPECTRUM OF LOW-MASS DRELL-YAN PAIRS 

In Drell-Yan production, if both, the physically measured dilepton mass M and the 
transverse momentum p T are large, the cross section in a collision between hadrons (or 
nuclei) A and B, A(P A ) + B{P B ) — > 7*(^ 11) + X, can be factorized systematically in QCD 
perturbation theory and expressed as [2] 

^APdydS = (SO 5 / dXl ^ A{Xl ^ ] I dx ^l^ ") ^J^f ^ X2 ' M ' PT > m ^ 

(1) 

The sum J2 a ,b runs over all parton flavors; (f> a /A and <p b / B are normal parton distributions; and 
\i represents the renormalization and factorization scales. The do ' a b->7* x / 1 dp\dy in Eq. (1) is 
the short- distance probability for partons of flavors a and b to produce a virtual photon of 
invariant mass M and is calculable perturbatively in terms of a power series in a s (/x). The 
scale /i is of the order of the energy exchange in the reaction, /i ~ y 1 M 2 + p\. 

The transverse momentum (jpx) distribution of the dileptons can be divided into three 
regions: low px (<^C M), intermediate pt (~ M), and high pr (>> M) regions. When both 
the physically measured M and pr are large and are of the same order, the short-distance 
partonic part do ' a b->7* x / 1 dp\dy in Eq. (1) can be calculated reliably in conventional fixed- 
order QCD perturbation theory in terms of a power series in a s (fj). However, when pt 
is very different from M, the calculation of Drell-Yan production in both low and high 
Pt regions becomes a two-scale problem in perturbative QCD, and the calculated partonic 
parts include potentially large logarithmic terms proportional to a power of ln(M/pr)- As 
a result, the higher order corrections in powers of a s are not necessary small. The ratio 
o NLO /o LO [oc a s x (large logarithms)] can be of order 1, and convergence of the conventional 
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perturbative expansion in powers of a s is possibly impaired. 

In the low p T region, there are two powers of ln(M 2 /p T ) for each additional power of 
a s , and the Drell-Yan pr distribution calculated in fixed-order QCD perturbation theory is 
known not to be reliable. Only after all-order resummation of the large a™ ln 2n+1 (M 2 /p^) 
do predictions for the p T distributions become consistent with data [3,4]. We demonstrate 
in Sec. Ill that low-mass Drell-Yan production at pr as low as Aqcd at LHC energies can 
be calculated reliably in perturbative QCD with all order resummation. 

When p T > M/2, the lowest-order virtual photon "Compton" subprocess: g + q — > 
7* + q dominates the px distribution, and the high-order contributions including all-order 
resummation of a" ln n_1 (M 2 /p^) preserve the fact that the p T distributions of low-mass 
Drell-Yan pairs are dominated by gluon initiated partonic subprocesses [5]. We show in 
Sec. IV that the px distribution of low-mass Drell-Yan pairs can be a good probe of the 
gluon distribution and its nuclear dependence. We give our conclusions in Sec. V. 

III. LOW TRANSVERSE MOMENTUM REGION 

Resummation of large logarithmic terms at low p T can be carried out in either p T or 
impact parameter (6) space, which is the Fourier conjugate of px space. All else being 
equal, the b space approach has the advantage that transverse momentum conservation is 
explicit. Using renormalization group techniques, Collins, Soper, and Sterman (CSS) [6] 
devised a b space resummation formalism that resums all logarithmic terms as singular as 
(1/Pt) \n m (M 2 /p T ) when p T — > 0. This formalism has been used widely for computations 
of the transverse momentum distributions of vector bosons in hadron reactions [7]. 

At low-mass M, Drell-Yan transverse momentum distributions calculated in the CSS b- 
space resummation formalism strongly depend on the non-perturbative parameters at fixed 
target energies. However, it was pointed out recently that the predictive power of perturba- 
tive QCD (pQCD) resummation improves with total center-of-mass energy (y/s), and when 
the energy is high enough, pQCD should have good predictive power even for low-mass 



3 



Drell-Yan production [3]. The LHC will provide us a chance to study low-mass Drell-Yan 
production at unprecedented energies. 

In the CSS resummation formalism, the differential cross section for Drell-Yan production 
in Eq. (1) is reorganized as the sum 

, j (resum) , (Y) 

d& AB->ll(M)X _ a<7 AB^ll(M)X a ° AB^ll(M)X 

dM 2 dydp 2 T ~ dM 2 dydp 2 T dM 2 dydp\ ' U 
The all-orders resummed term is a Fourier transform from the 6-space, 

7 (resum) 

du 1 r - 1 r 

dijtfgf = p^/^^'H-'IM^,,^) = -JdbMrrVbW^M,^) , (3) 



where J is a Bessel function, xa = e y M/y/s and = e~ y M/y/s, with rapidity y and 
collision energy y/s. In Eq. (2), the <j( rcsum ) term dominates the px distributions when 
Pt <C M, and the c/ y ) term gives corrections that are negligible for small p?, but become 
important when px ~ M. 

The function VT(6, M, xa, %b) resums to all orders in QCD perturbation theory the sin- 
gular terms that would otherwise behave as S 2 (px) and (l/p?) ln m (M 2 /p^) in transverse 
momentum space, for all m > 0, and can be calculated perturbatively for small b , 

W(b, M, x A , x B ) = e~ s{b ' M) W(b, c/b, x A , x B ) = W peTt (b, M, x A , x B ) , (4) 

where all large logarithms from ln(c 2 /6 2 ) to ln(M 2 ) have been completely resummed into the 
exponential factor S(b,M) = J^ b2 dfi 2 / fi 2 [ln(M 2 / fi 2 )A(a s (fi)) + B(a s (fi))] with functions 
A and B given in Ref. [6], and c = 2e _TB with Euler's constant 7# « 0.577. The function 
W(b,c/b,XA,XB) in Eq. (4) is given in terms of modified parton distributions from hadron 
A and B [6]. With only one large momentum scale 1/6, W(b,c/b,XA,xs) is perturbatively 
calculable. Since the perturbatively resummed W pcrt (b, M, xa, xb) in Eq- (4) is only reliable 
for the small b region, an extrapolation to the nonperturbative large b region is necessary in 
order to complete the Fourier transform in Eq. (3). 

In the original CSS formalism, a variable 6* and a nonperturbative function 
Fcg S {b, M,xa,x b ) were introduced to extrapolate the perturbatively calculated W pert into 
the large b region such that the full 6-space distribution was of the form 



W css (b, M, x A , x B ) = W pcrt (K, M, x A , x B ) F^ s p s (b, M, x A , x B ) , (5) 

where 6* = b/yjl + (b/b max ) 2 , with b max = 0.5 GeV -1 . This construction ensures that 
b* < b max for all values of b. 

In terms of the b* formalism, a number of functional forms for the F^ p s have been 
proposed. A simple Gaussian form in b was first proposed by Davies, Webber, and Stirling 
(DWS) [8], 

F Dws(b, M,x A ,x B ) = exp {-(<?! +g 2 \n(M/2M ))b 2 }, (6) 

with the parameters M = 2 GeV, g-y = 0.15 GeV 2 , and g 2 = OA GeV 2 . In order to take into 
account the appearent dependence on collision energies, Ladinsky and Yuan (LY) introduced 
a new functional form [9], 

Fg?(b, M, x A , x B ) = exp {-((?!+ g 2 ln(M/2M ))6 2 - 9l g 3 \n(100x A x B )b} , (7) 

with M = 1.6 GeV, g x = O.lllg;^ GeV 2 , g 2 = 0.58^ GeV 2 , and ^ 3 = — 1.5±g;l GeV" 1 . 
Recently, Landry, Brook, Nadolsky, and Yuan proposed a modified Gaussian form [10], 

Fblny 

(6, M, x A , i B ) = exp {- [ 9l + ^ 2 ln(M/2M ) + gi g 3 ln(100 MB )] 6 2 } , (8) 

with M = 1.6 GeV, 9l = 0.21+°;^ GeV 2 , g 2 = 0.68±g;^ GeV 2 , and # 3 = -0.6^. All 
these parameters were obtained by fitting low energy Drell-Yan and high energy W and Z 
data. Note, however that the 6* formalism introduces a modification to the perturbative 
calculation, and the size of the modifications strongly depends on the non-perturbative 
parameters in F NP (b, M,xa,x b ), the intermediate boson mass M, and collision energy ^Js 
[7]- 

A remarkable feature of the 6-space resummation formalism is that the resummed expo- 
nential factor exp[— S(b, M)] suppresses the 6-integral when b is larger than 1/M. It can be 
shown using the saddle point method that, for a large enough M, QCD perturbation theory 
is valid even at pt = [6]. For high energy heavy boson (W, Z, and Higgs) production, the 
integrand of 6-integration in Eq. (3) at pr = is proportional to bW(b, Q, x A , x B ), which 
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has a saddle point b sp well within the perturbative region (b sp < b max ), and therefore, the 
6-integration in Eq. (3) is dominated by the perturbatively resummed calculation. The un- 
certainties from the large-6 region have very little effect on the calculated pr distributions, 
and the resummation formalism is of a good predictive power. 

On the other hand, in the low energy Drell-Yan production, there is no saddle point in the 
perturbative region for the integrand in Eq. (3), and therefore the dependence of the final 
result on the non-perturbative input is strong [3,4]. However, as discussed in Ref.s [3,7,11], 
the value of the saddle point strongly depends on the collision energy y/s, in addition to its 
well-known M 2 dependence. 

Figure 1 shows the integrand of the ^-integration in Eq. (3) at pr = for production 
of Drell-Yan pairs of mass M = 5 GeV in proton-proton collisions at y^s = 5.5 TeV and 
bmax = 0.5 GeV -1 . Different curves represent different extrapolations to the large-6 region. 
The three curves (dashed, dotted, and dot-dashed) are evaluated using the 6* formalism 
with the DWS, LY, and BLNY nonperturbative functions, respectively. Although these 
three nonperturbative functions give similar 6-space distributions for heavy boson produc- 
tion at Tevatron energies, they predict very different 6-space distributions for low-mass 
Drell-Yan production at LHC energies even within the perturbative small-6 region. Since 
the 6-distribution in Fig. 1 completely determines the resummed pt distribution through 
the b- integration weighted by the Bessel function J (p T b), we need to be concerned with the 
uncertainties of the resummed low-mass px distributions calculated with different nonper- 
turbative functions. We use the CTEQ5M parton distribution function [12] throughout. 

In order to improve the situation, a new formalism of extrapolation (QZ) was proposed 

[3], 



W(b, M,x A ,x B ) 



WP ert (b,M,x A: x B ) b<b max 

> ( 9 ) 

W peTt (b max , M, x A , x B ) F NP (b, M, x A , x B ; b max , a) b> b max 



where the nonperturbative function F NP is given by 

F NP = exp {- HM%lJc 2 ) [ 9l ((b 2 r - (b 2 max ) a ) +92 (b 2 - b 2 max )} - g 2 (b 2 - b 2 nax ) } . (10) 



Here, b max is a parameter to separate the perturbatively calculated part from the non- 
perturbative input, and its role is similar to the b max in the 6* formalism. The term pro- 
portional to g\ in Eq. (10) represents a direct extrapolation of the resummed leading power 
contribution to the large b region; and the parameters g\ and a are determined by the con- 
tinuity of the function W (6, M, xa, %b) at b max . On the other hand, the values of g 2 and g 2 
represent the size of nonperturbative power corrections. Therefore, sensitivity on the g 2 and 
g 2 in this formalism clearly indicates the precision of the calculated px distributions. 

The solid line in Fig. 1 is the result of the QZ parameterization with b max = 0.5 GeV -1 
and g 2 = g 2 = 0. Unlike in the 6* formalism, the solid line represents the full perturbative 
calculation and is independent of the nonperturbative parameters for b < b max . The differ- 
ence between the solid line and the other curves in the small b region, which can be as large 
as 40%, indicates the uncertainties introduced by the 6* formalism. 

It is clear from the solid line in Fig. 1 that there is a saddle point in the perturbative 
region even for the dilepton mass as low as M = 5 GeV in Drell-Yan production at y/s = 5.5 
TeV. At that energy, xa, %b ~ 0.0045. For such small values of x, the PDFs have very strong 
scaling violation, which leads to a large parton shower. It is the large parton shower at the 
small x that strongly suppresses the function W(b, c/b, xa, xb) in Eq- (4) as b increases. 
Therefore, for a b max ~(a few GeV) -1 , the predictive power of the 6-space resummation 
formalism depends on the relative size of contributions from the small-6 (b < b max ) and 
large-6 {b > b rnax ) regions of the 6- integration in Eq. (3). With a narrow b distribution 
peaked within the perturbative region for the integrand, the 6-integration in Eq. (4) is 
dominated by the small-6 region, and therefore, we expect pQCD to have good predictive 
power even for low M Drell-Yan production at LHC energies. 

Figure 2 presents our prediction to the fully differential cross section da/dM 2 dydp^ for 
Drell-Yan production in pp collisions at y/s = 5.5 TeV. Three curves represent the different 
order of contributions in a s to the perturbatively calculated functions A(a s ), B(a s ), and 
C(a s ) in the resummation formalism [13]. The solid line represents a next-to-next-to-leading- 
logarithmic (NNLL) accuracy corresponding to keeping the functions, A(a s ), B(a s ), and 
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C(a s ) to the order of oP s , a 2 , and a], respectively. The dashed line has a next-to-leading- 
logarithmic (NLL) accuracy with the functions, A(a s ), B(a s ), and C(a s ) at the a 2 , a s , and 
a", respectively, while the dotted line has the lowest leading-logarithmic (LL) accuracy with 
the functions, A(a s ), B(a s ), and C(a s ) at the a s , and ct°, respectively. Similar to what 
was seen in the fixed order calculation, the resummed pt ditribution has a K-factor about 
1.4-1.6 around the peak due to the inclusion of the coefficient 

In Eq. (10), in addition to the g\ term from the leading power contribution of soft 
gluon showers, the g 2 term corresponds to the first power correction from soft gluon showers 
and the g 2 term is from the intrinsic transverse momentum of the incident parton. The 
numerical values of g 2 and g 2 have to be obtained by fitting the data. From the fitting 
of low energy Drell-Yan data and heavy gauge boson data at the Tevatron, we found that 
the intrinsic transverse momentum term dominates the power corrections and it has a weak 
energy dependence. For convenience, we combine the parameters of the b 2 term as G 2 = 
g 2 \n(M 2 b 2 nax /c 2 ) + g 2 . For M = 5 GeV and y = 0, we use G 2 ~ 0.25 in the discussion here 
[3]. To test the G 2 dependence of our calculation, we define 



where the numerator represents the result with finite G 2 , and the denominator contains no 
power corrections (G 2 = 0). 

The result for R G2 is shown in Fig. 3. R G2 deviates from unity less than 1%. The 
dependence of our result on the non-perturbative input is indeed very weak. 

Since the G 2 terms represent the power corrections from soft gluon showers and partons' 
intrinsic transverse momentum, the smallness of the deviation of Rg 2 from unity also means 
that leading power contributions from gluon showers dominate the dynamics of low-mass 
Drell-Yan production at LHC energies. Even though the power corrections will be enhanced 
in nuclear collisions, we expect it to be still less than several percent [11]. The isospin effects 
are also small here, because xa and xb are small. 

Since the leading power contributions from initial-state parton showers dominate the 




(11) 
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production dynamics, the important nuclear effect is the modification of parton distribu- 
tions. Because the xa and x B are small for low- mass Drell-Yan production at LHC energies, 
shadowing is the only dominant nuclear effect. In order to study the shadowing effects, we 
define [11] 

„ / _ da ( ^ lT(M)x (p T ,Z A /A,Z B /B) I da AB ^ lI(M)x (p T ) 

Ksh{ P T >~ dM^dydp 2 T I dM*dydp 2 T ' 1 ] 

We plot in Fig. 4 the ratio R s h as a function of pr in pPb and PbPb collisions at y/s = 5.5 TeV 
for M = 5 GeV and y — 0. The EKS parameterizations of nuclear parton distributions [14] 
were used to evaluate the cross sections in Eq. (12). Fig. 4 shows that R s h decreases about 
30% from pPb to PbPb collisions. It is clear that low-mass Drell-Yan production at pPb 
and PbPb can be a good probe of nuclear shadowing. 

IV. HIGH TRANSVERSE MOMENTUM REGION 

The gluon distribution plays a central role in calculating many important signatures 
at hadron colliders because of the dominance of gluon initiated subprocesses. A precise 
knowledge of the gluon distribution as well as its nuclear dependence is absolutely vital for 
understanding both hard and semihard probes at LHC energies. 

It was pointed out recently that the transverse momentum distribution of massive lepton 
pairs produced in hadronic collisions is an advantageous source of constraints on the gluon 
distribution [15], free from the experimental and theoretical complications of photon isolation 
that beset studies of prompt photon production [16,17]. Other than the difference between 
a virtual and a real photon, the Drell-Yan process and prompt photon production share the 
same partonic subprocesses. Similar to prompt photon production, the lowest-order virtual 
photon "Compton" subprocess: g+q — > 7* +q dominates the Pt distribution when p T > M/2, 
and the next-to-leading order contributions preserve the fact that the px distributions are 
dominated by gluon initiated partonic subprocesses [15]. 

There is a phase space penalty associated with the finite mass of the virtual photon, and 
the Drell-Yan factor a em /(37rM 2 ) < 1CT 3 /M 2 in Eq. (1) renders the production rates for 
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massive lepton pairs small at large values of M and pr- In order to enhance the Drell-Yan 
cross section while keeping the dominance of the gluon initiated subprocesses, it is useful 
to study lepton pairs with low invariant mass and relatively large transverse momentum 
[5]. With the large transverse momentum p T setting the hard scale of the collision, the 
invariant mass of the virtual photon M can be small, as long as the process can be identified 
experimentally, and the numerical value M ^> Aqcd- For example, the cross section for 
Drell-Yan production was measured by the CERN UA1 Collaboration [18] for virtual photon 
mass M e [2m M , 2.5] GeV. 

When p\ M 2 , the perturbatively calculated short-distance partonic parts, 

da a b-^ 1 *x/dp%dy in Eq. (1), receive one power of the logarithm \n(p^/M 2 ) at every order 
of a s beyond the leading order. At sufficiently large p Tl the coefficients of the perturbative 
expansion in a s will have large logarithmic terms, and these high order corrections may not 
be small. In order to derive reliable QCD predictions, resummation of the logarithmic terms 
\n m (p^/M 2 ) must be considered. It was recently shown [5] that the large ln m (p^/M 2 ) terms 
in low-mass Drell-Yan cross sections can be systematically resummed into a set of pertur- 
batively calculable virtual photon fragmentation functions [19], and similar to Eq. (2), the 
differential cross section for low-mass Drell-Yan production at large px can be reorganized 
as 

aa AB^U(M)X _ uu AB^ll(M)X uu AB^ll(M)X ,..„,. 

dM 2 dydp 2 T ~ dM 2 dydp 2 T dM 2 dydp 2 T ' [ } 

where <j( rcsum ) includes the large logarithms and can be factorized as [5] 

^^_(rcsum) 

dWdydft = (SO J c fteih/Aixuri j dx^/six^) 

x / ^ ^§(pc t = pt/z) /Wx(^,/4; m 2 ) , (14) 



with the factorization scale p, and fragmentation scale /xp, and the virtual photon fragmenta- 
tion functions D c ^ 7 * (z, p 2 F ; Q 2 ). The <j( Dir ) term plays the same role as term in Eq. (2), 
and it dominates the cross section when p T — > M. 
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Figure 5 presents the fully resummed transverse momentum spectra of low-mass Drell- 
Yan production in pp collisions with M = 5 GeV at y = and y/s = 5.5 TeV (solid). 
For comparison, we also plotted the leading order spectra calculated in conventional fixed 
order pQCD. The fully resummed distribution is larger in the large pr region and smoothly 
convergent as pr — > 0. In addition, as discussed in Ref. [5], the resummed differential 
cross section is much less sensitive to the changes of renormalization, factorization, and 
fragmentation scales, and should be more reliable than the fixed order calculations. 

To demonstrate the relative size of gluon initiated contributions, we define the ratio 

_ doAB^r (M)x (gluon-initiated) / da A B^ 1 *{M)x 
9 dp\dy I dp\dy 

The numerator includes the contributions from all partonic subprocesses with at least one 
initial-state gluon, and the denominator includes all subprocesses. 

In Fig. 6, we show R g as a function of pt in pp collisions at y = and = 5.5 TeV 
with M = 5 GeV. It is clear from Fig. 6 that gluon initiated subprocesses dominate the low- 
mass Drell-Yan cross section and that low-mass Drell-Yan lepton-pair production at large 
transverse momentum is an excellent source of information on the gluon distribution [5]. 
The slow falloff of R g at large p T is related to the reduction of phase space and the fact that 
cross sections are evaluated at larger values of the partons' momentum fractions. 



V. CONCLUSIONS 

In summary, we present the fully differential cross section of low-mass Drell-Yan pro- 
duction calculated in QCD perturbation theory with all-order resummation. For p T <C M, 
we use CSS 6-space resummation formalism to resum the large logarithmic contributions as 
singular as ln m (M 2 /Pt)/Pt t° a ^ orders in a s . We show that the resummed p T distribution 
of low-mass Drell-Yan pairs at LHC energies is dominated by the perturbatively calculable 
small 6-region and thus reliable for pt as small as Aqcd- Because of the dominance of small 
x PDFs, the low-mass Drell-Yan cross section is a good probe of the nuclear dependence 
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of parton distributions. For px ^ M, we use a newly derived QCD factorization formal- 
ism [5] to resum all orders of \n m (p T /M 2 ) type logarithms. We show that almost 90% of 
the low-mass Drell-Yan cross sections at LHC energies is from gluon initiated partonic sub- 
processes. Therefore, the low-mass Drell-Yan cross section at pr > M is an advantageous 
source of information on the gluon distribution and its nuclear dependence — shadowing. 
Unlike other probes of gluon distributions, low-mass Drell-Yan does not have the problem 
of isolation cuts associated with direct photon production at collider energies, and does not 
have the hadronization uncertainties of J /if; and charm production. Moreover, the precise 
information on dilepton production from the Drell-Yan channel is critical for studying charm 
production at LHC energies. 
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FIG. 1. The 6-space resummed functions bW(b) in Eq. (3) for Drell-Yan production of dilepton 
mass M = 5 GeV at ^/s = 5.5 TeV with the ZQ (solid), DWS (dashed), LY (dotted), and BNLY 
(dotdashed) formalism of nonperturbative extrapolation. 
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FIG. 2. Differential cross section da / dM 2 dydp\ for production of Drell-Yan pairs of M = 5 GeV 
in pp collisions at the LHC with y = 0, and yfs = 5.5 TeV with the NNLL (solid), NLL (dashed), 
and LL (dotted) accuracy. 



14 




Pt (GeV) 

FIG. 3. The ratio R G2 defined in Eq. (11) with G 2 = 0.25 GeV 2 for the differential cross section 
shown in Fig. 2. 
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FIG. 4. R s h as a function of pr using EKS shadowing in PbPb collisions at y/s = 5.5 TeV. 
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FIG. 5. Differential cross section da I dM 2 dydp\ for production of Drell-Yan pairs of M = 5 GeV 
in pp collisions at y/s = 5.5 TeV with low and high px resummation (solid), in comparison to 
conventional lowest result (dashed). 
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FIG. 6. Ratio of gluonic over total contributions to Drell-Yan production at the LHC, R g , 
defined in Eq. (15) with M = 5 GeV at ^ = 5.5 TeV. 
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